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MODE I ANALYSIS OF A CRACKED CIRCULAR DISK SUBJECT 


TO A COUPLE AND A FORCE 
by Bernard Gross 
Lewis Research Center 


SUMMARY 

Mode I stress intensity coefficients were obtained for an edge- 
cracked disk (Round Compact Specimen) . Results for this plane elasto- 
static problem, obtained by a boundary collocation analysis are presented 
for ratios 0.35 < A/D < 1, where A is the crack length and D is the 
disk diameter. The results presented are for two complementary types of 
loading. By superposition of these results the stress intensity factor 
Kj., for any practical load line location of a pin- loaded round compact 

specimen can be obtained. 


INTRODUCTION 


The round compact specimen <edge-cracked disk) is currently being 
considered by ASTM Committee E-24 on Fracture Testing for incorporation 
into ASTM Standard Method of Test E 399 on Plane Strain Fracture Tough- 
ness of Metallic Materials. In reference 1 values of the mode I stress 
intensity coefficients Y and r M for round compact specimens are 


given for the range of ratio of crack length to specimen diameter, A/D, 
from 0.65 to 0.9 (where crack length A here is the distance from the 
crack tip to the specimen circumference). The coefficients and I\ r 

apply to two complementary types of specimen loading; T to a nominal 


uniform distribution of stress across the net section, and to a 


nominal bending stress distribution (fig. 1). While these types of load- 
ing are impractical in themselves, the two coefficients can be combined 
to represent any practical case of loading of the specimen by a pair of 
equal and opposite forces normal to the crack (pin loading). The appro- 
priate value of the stress intensity factor is then obtained from this 
combination of coefficients. The approach used in reference 1 was to 
obtain values of V and by the boundary collocation method for a 

set of ring segment specimens (ref. 2) with successively diminished ratios 
of internal to external radius to a minimum ratio of 0.07. These results 
were then extrapolated to estimate values for the limit case of zero 
internal radius, corresponding to the edge-cracked disk. 
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To extend and assess the results of reference l t a more direct 
approach has now been taken which employs boundary value stress functions 
and derivatives appropriate to the cracked disk. The present results 
cover the range of ratio of crack length to specimen diameter, A/D, from 
0.35 to 1. These results are given both in tabular form and in the form 
of interpolation equations obtained by multiple linear and nonlinear re- 
gression analysis. It should be appreciated that there is a limitation 
on the applicability of the present results to practical pin-loaded 
cracked disk specimens. When the crack tip is close to the load line 
there will be a significant difference between the actual distribution 
of loading forces and that assumed in the present model, and the values 
of the corresponding stress intensity factors will be significantly dif- 
ferent. It is estimated that this effect will be negligible when the 
distance of the crack tip from the load line exceeds 0.15 D. 


LIST OF SYMBOLS 


a 

B 

D 

K I 

L 

M 

P 

W«(D/2+L) 


crack length measured from crack tip to the circum- 
ference 

crack length measured from the crack tip to load line 
specimen thickness 
disk diameter 

Mode I stress intensity factor K T - K TW when a - 0, 

1 LM p 

and K_ » K T when o x . * 0 
I Ip M 

load line location 

resultant moment at nominal neutral axis of specimen 
applied pin load 

distance measured from load line to circumference of 
specimen 


F 7o p -hJ M ) VA(1-A/D) 8treSS inten8ity “efficient r - T p when c M 
a M -6M/[B(D-A) 2 ] 


0 , 


and r - when a * 0 
M p 


a -P/[B(D-A)] 


component of fictitious normal net stress at the crack 
tip due to M 

component of fictitious normal net stress due to load P 
stress function 
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APPROACH 


As shown in figure 1, the cracked round specimen has a diameter 
and is loaded through pins by a couple Mo and opposed tensile forces 
P normal to the crack with load line offset a distance L from the 
center of the specimen. Calculations were performed for each of two 
complementary types of loading: case 1, nominal uniform tension and 

case 2, nominal bending. For case 1, provided the ratio a/D (“A/D 
+L/D-1/2) is at least 0.15, and regardless of the load line location 
the computed stress intensity factor is a function of two vari- 


D 


L 


» 


ables: nominal uniform stress and crack length ratio A/D. This 

case is obtained by adjusting the boundary load condition along @ ® 
to produce at the nominal neutral axis position (D-A)/2, a resultant 
moment of zero, and force equal to the applied load P. 


For case 2, provided the ratio a/D is at least 0.15, and regard- 
less of the load line location L the computed stress intensity factor 
is a function of two variables: nominal bending stress at the 

crack tip a^, and crack length ratio A/D. This case is obtained by 

applying only a cour le at boundary 


Stress intensity coefficienty 




and 


f are defined as 


r 

p 




sfia - a / d ) 







0 M ^A(l - A/D) 




( 1 ) 


r 


(a p + c M ) tykti - A/D)" J 


For the general case (fig. 1) , there are two independent variables L/D 
and A/D. As an illustration, for a disk with load P and couple M q 
at load line location L: 


P 

°p “ B(D - A) 


and 
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3£2M o + P(2L + A)] 
B(D - A) 2 


By superposition we have 


K I “ K IP + k im 


( 2 ) 


Applying equation (1) and equation (2) we obtain 


W rp * (', " •«) 


M 


(3) 


Replacing a and o w by their equivalent definitions we obtain 

p M 


P(D - A)T p + 3[2M o + P(2L + A) ]r^ 
P(D - A) + 3[2M + P(2L + A)] 


(4) 


If at loadline L, a load P is applied through pins and the moment 
M ■ 0, after algebraic manipulation of equation (4) we obtain a more 

common form 


KB V D 
P 


4 


A/D 


1 - (A/D) 


r + 3 (A/D) + (2L/D) " 

P 1 - (A/D) ‘M 


(5) 


RESULTS AND DISCUSSION 

As shown in figure 1, the resultant solution of the cracked disk 
problem is obtained by combining two complementary types of loading. 

The first type of loading is based on a nominal constant mid-net 
section stress * P/[B(D-A)] resulting in the stress intensity coef- 
ficient T as a function of a and A/D. The second solution is 
P P 

based on a nominal pure bending stress at the mid net section where 
a M*3^ (A+2L)+M ] / [B(D-A) from which the stress intensity coefficient 

r u as a function of cl. and A/D is obtained. 

T and IY, are obtained using the boundary collocation method 
p M. 

with 60 boundary stations and an over determined system of equations as 
detailed in reference 3. Twenty equally spaced boundary stations were 
taken along the straight portion of the boundary @ ® and 40 stations 
equally spaced were taken along the curved boundary nn © (fig. Al). 
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The appropriate stress function boundary conditions for both cases are 
given in the appendix. The limit values for and V M as A/D + l were 

obtained from reference 4. The values of r and given in table I 

p M 

were obtained for a load line offset L - 0.35D. The values given in 
table II of T were obtained directly by boundary collocation with the 
load line L * D/4. 

It has been shown (ref. 5) that the pin loaded holes can have a 
significant effect on the stress intensity factor. It is therefore 
estimated that the results presented herein will apply when the distance 
from the crack tip to the load line exceeds 0.15 D. 

Table II (eq. (4) with M q - 0) contains a comparison of the present 

results with those of references 1, 6, and 7 , for a pin loaded cracked 
disk with load line at L * D/4. Excellent agreement is obtained. In- 
cluded in this table are the stress intensity coefficients of the stand- 
ard rectangular compact specimen (ref. 8). A geometric comparison be- 
tween the standard rectangular and round compact specimen is given in 
figure 2. Fitting functions for and were obtained by linear 

and nonlinear least squares best fit regression analyses respectively. 
These functions are: 

r M - 5.611 - 21.246 £ + 38.149 - 12-169 (j) + 10.320 

and 


r 

p 


0.501 + 


0.02104 [(A/D) - 0.5251 

[(A/D) - 0.525] [(A/D) - 0.0798] +0.0625 


These functions are considered to be accurate to less than a percent of 
the computed solution in the range 0.35 < A/D < 1. Clearly, the accuracy 
of application is dependent upon how well the assumed model boundary 
conditions approximate the real boundary conditions. 


EXAMPLE 

Two examples now follow which will demonstrate the use of table I 
in conjunction with equation (4) and equation (5). 


Example 1 

For a pin loaded disk with load P at load line location L - D/4, 
to find the value of KB h! W/P where W ■ 0.75 D and a/W ■ 0.4. Thus 


6 


a - A - 0,25 D, A/D - 0.55 and from table I, r p ■ 0.509 and * 1.064. 
From equation (5) since M - 0, we have KBVWp - V 574 KBV57p - 7.618. 

There Is a small difference between this value and the value 7.613 
given In table XX. The values In table XX are obtained by direct appli- 
cation of the boundary conditions for this loading and the value 7.618 
In the above example Is obtained by combination. 


Example 2 

Given a disk with crack length A, load P 9 and couple M q at 

L ■ 0 9 to find the resultant stress intensity coefficient T. From 

table X we obtain T and and using equation (4) we have 

p m 

P(D - A)r + 3(2M + ?A)r M 
p O n 

PD + 6M + 2PA 
o 


r 
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APPENDIX 

The results presented here were obtained by boundary collocation 
analysis of a homogeneous isotropic specimen under plane elastostatlc 
conditions. This method is described in detail in reference 3. The 
necessary boundary conditions to be satisfied by the stress function 
and its normal derivative (fig* Al) are as follows 

Along ® (D 


pP(j> 


/ -sin 
\2a + 


M 

a sin 4 ■ cos a cos b \ . o/sin 2<fr - 2$ cos 2a _ 

** - cos 2a/ 2 \2a cos 2a - sin 2a ~ / 


sin 2a 


2a 


+ pp cos $ 


- Pp sin $ 


3 

-sin a 


sin 2a 


\2a + 

( a sin a + 
2a + 


a cos a - sin a cos 


2a - sin 2a 


cos a sin a 
sin 2a 


cos 


2a - sin 
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a \ 

n 2a J 


and 


eos (♦ + a) - sin (* + a) 


Along ®© 


X ■ P ^ (cos ♦ + sin 2a) + Pe - M o 


and 


- P cos $ 

The symbols p t <f>, a, and e are defined in figure Al. 
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TABLE I. - VALUES OF STRESS INTENSITY COEFFICIENTS 
T p AND r M FOR THE EDGE CRACKED DISK 


WHEN a/D > 0.15 


1 

A/D 

r F 

r M 

0.35 

0.2587 

1.6157 

.40 

.3857 

1.4172 

.45 

.4538 

1.2695 

.50 

.4901 

1.1551 

.55 

.5090 

1.0637 

.60 

.5184 

.9865 

.65 

.5228 

.9255 

.70 

.5246 

.8717 

.75 

.5251 

.8252 

.80 

.5251 

.7845 

.85 

.5249 

.7487 

.90 

.5242 

.7168 

.95 

.522 

.682 

*1.00 

.521 

.663 


a 

Limit values obtained from ref. 4. 
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TABLE II. - COMPARISON OF PRESENT STRESS INTENSITY PARAMETER 
KjBVw/P AS A FUNCTION OF a/W WITH THOSE RESULTS OF 
REFERENCES 1 # 6, 7, AND 8 (AS SHOWN IN FIG. 2 


FOR L/D - 0.25 AND W - 0.75 D) 


0.200 

.267 

,300 

.333 

.400 

.467 

,500 

.533 

.600 

.667 

.700 

.733 

.800 

.867 

.933 


Direct 

Collocation 

Experi- 

Finite 

Compliance 

Standard 

collo- 

ring seg- 

mental 

elements 

resvlts 

rectangular 

cation 

ment 

results 

Ref. 7 

Ref. 7 

compact 

solid 

Ref. 1 

Ref. 6 



specimen 

disk 





Ref. 8 

K BA/W/P 


4.745 

5.472 





6.365 

7.613 

9.203 

— — — — — 

6.244 

7.299 

8.948 

11.378 

14.532 

19.368 

11.926 

14.820 

19.518 

11.314 

14.721 

19.680 

27.411 

42.773 

79.542 

226. 

27.490 
42. 745 
79.523 

26.879 


— 


— 

— 

4.750 



5.844 



— 

— 

7,333 

.10,185 

10.286 

9.631 

14.738 

14.566 

13.62 

22.653 

22.798 

21.56 

42.542 


41,07 




















Aj 

(wm ft# 




finite fat# sfie/c. Conf**'to*J or r&£ A(MO CtwAHr 

#tr& smv&he* /eecnwf***** ctwAKr *>* 

i*r» ts/r*. feAfrv*/ L * Dfo 


